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A  FORTRAN  COMPUTER  PROGRAM  TO  COMPUTE 
THE  RADIATION  PATTERN  OF  AN  ARRAY-FED 
PARABOLOID  REFLECTOR 


INTRODUCTION 

Computing  the  radiation  pattern  of  a  paraboloid  reflector  has  been  treated  widely  in  literature  [1- 
4];  however,  most  of  this  treatment  concerns  only  the  problem  of  single-feed  source  whether  the 
source  is  at  the  focal  point  or  offset.  In  all  these  treatments,  when  the  source  is  offset,  it  is  always 
assumed  that  the  offset  is  very  small  and  that  the  distance  from  the  source  to  a  point  on  the  reflector 
surface  can  be  approximated  by  adding  a  few  linear  terms  to  the  case  of  the  feed  at  the  focal  point. 
This  approximation  may  be  acceptable  in  a  single-source  case  when  the  feed  point  usually  is  located 
close  to  the  focal  point  to  avoid  high  sidelobes.  However,  in  the  case  of  an  array-fed  system,  this 
assumption  may  not  be  valid. 

To  facilitate  fast  computation,  the  longitudinal  current  component  is  ignored  in  this  computation 
[1]  The  error  introduced  by  this  omission,  in  general,  is  small  in  the  region  close  to  the  main  beam. 
However,  when  fields  far  from  the  main  beam  are  of  interest,  the  amount  of  phase  error  can  not  be 
ignored.  Galindo-Israel  and  Mittra  [31  introduced  a  method  to  correct  this  error;  however,  this  correc¬ 
tion  depends  on  the  location  of  the  offset  feed.  It  probabily  would  not  be  efficient  for  a  multiple-feed 
case. 


In  this  report  a  computer  program  to  compute  the  radiation  pattern  of  an  array-fed  paraboloid 
reflector  is  presented.  This  program  is  written  so  that  it  is  efficient  for  multiple-feed  sources  and  the 
solution  will  be  exact.  No  simplification  or  other  approximating  assumption  is  made. 

The  purpose  of  developing  this  computer  program  is  threefold. 

•  It  will  be  used  as  a  reference  to  check  the  validity  and  efficiency  of  other  programs  which 
may  be  developed  using  a  more  simplified  approach. 

•  It  is  an  interim  step  to  develop  a  computer  program  to  compute  the  radiation  pattern  of  a 
dual-reflector  system  which  is  currently  under  a  feasibility  study  at  NRL. 

•  This  program  may  be  used  to  investigate  the  pattern-synthesis  method  for  an  array-fed 
reflector  antenna  system. 

MATHEMATICAL  FORMULATION 

Figure  1  shows  an  array-fed  paraboloidal  reflector  antenna  system.  The  radiation  field  at  a  point  $ 
and  9  has  the  form 

!;(♦,  9)  -  C  f  T  7„(e,  4)  -r  exp  l-jkiph  -  p  •  A)]  da,  (1) 

-  Pm 


Manuscript  approved  May  9,  19*3. 
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Fig.  I  —  Configuration  of  an  array-fed 
parabaloid  reflector 
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where  Jj(9,  <f>)  is  the  current  density  induced  on  the  surface  of  the  reflector  and  p'„  is  the  distance  for 
the  nth  array  feed  element  to  a  point  on  the  surface,  and 

p'n-p-  e*.  (2) 

where  e„  is  a  position  vector  from  the  reference  point  to  the  location  of  the  nth  array  element.  For 
convenience,  the  reference  point  is  assumed  to  be  at  the  focal  point  of  the  parabola.  Unit  vector  R 
points  the  direction  from  the  current  source  to  a  field  point  which  is  a  function  of  <t>  and  9. 

Assume  that  only  part  of  the  paraboloid  is  used.  This  used  part  is  known  to  be  a  circle  at  the 
upper  half  of  the  reflector  (see  Fig.  1).  The  center  of  this  circle  is  at  a  distance  dt  above  the  r-axis. 
Any  point  on  this  circle  can  be  represented  by 

x  -  p  sin  0  cos  <t»  -  r  cos  <f>  (3a) 

y  -  p  sin  0  sin  ^  -  r  sin  4> ■  (3b) 

In  terms  of  the  coordinates  of  the  4>'  and  rt  of  this  circle, 

x  -  d\  +  cos  d>' 


and 


is 


(4) 

y  -  rj  sin 

The  position  vector  from  a  feed  element  in  the  array  to  a  point  (jr,  y,  z)  on  the  reflector  surface 


P* - TJi—  x  +  -  |_7p"  y  +  |p,[  * 


f  *  I  p  A  ,  I  A 

-  +  pvy  +  p«z- 

The  paraboloid  surface  can  be  described  by  the  following  equation: 


(5) 

(6) 


x2  +  y2  -  4f(f  +  z) 
where  /is  the  focal  length  of  the  parabola. 

The  surface  currents  7„  are  related  to  the  feed  element  exciting  field  by  the  following  relation: 

X  -  "  x  ?!  x  e„.  (7) 

where  e„  is  a  unit  vector  in  the  direction  of  the  tangential  field  on  the  surface  of  the  paraboloid  which  is 
excited  by  the  nth  array  element  and 


p'j\p;  X  X  j>;l 


(8) 


where  u„  is  a  unit  vector  in  the  direction  of  the  field  in  the  nth  array  element.  Assume  that  the  array 
elements  and  paraboloid  have  the  same  coordinates,  and  that 


«»  “  u„x  +  uvy  +  u„z 


(9) 


where  |«.|  —  1. 


The  normal  vector  n  in  Eq.  (7)  for  a  paraboloid  surface  has  the  following  form: 

2  j 


A  X  A  V  A  ,  A 

—  Ux--Q,  +  Ti[‘ 

—  +  n,  x  +  ity}  +  rtjZ. 


(10) 
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Aopendix  A  shows  that 

Jn  -  —7—'  ==-  U-  “n*(nyp'„y  +  n2p'm)  +  u^p'mny  +  u^rijx 
V  1  -  a2 

+  l UnxPnyn*  ~  “ny^xP'nx  +  "zP*)  + 


+  l"«P«"x  +  Uny"yPn2  ~  ««<«xP«  +  nyp'ny))z]. 


(11) 


where  a  -  fa  •  u„  -  p'm  +  unyp'ny  +  u^. 

The  u„,  and  um  are  the  normalized  excitation  field  at  the  nth  element  in  the  feed  array.  Usu¬ 
ally,  the  elements  are  linear  polarized  either  in  the  x  or  in  the  y  direction.  In  this  case  only  one  com¬ 
ponent  will  be  equal  to  unity. 


The  p'„  components  are 


P'nx  “  +  z,  cos  <*>'-  «„)/lp«l> 

(12a) 

J 

Pny  “  (ft  sin  4>'  ~  «v)/|p,;|. 

(12b) 

Pnz  "  I  - 

~  [4/2  -  (d,2  -  2d,zi  cos  <f>'  +  z,2)  1  -  e„  /  |pj. 

(12c) 

and  components  of  h  vector  are 

nx  -  -  (dx  +  zj  cos  ^')/U  I, 

(13a) 

ny  -  -  z,  sin <6 '/Inf, 

(13b) 

nt  -  2//UI, 

(13c) 

where  Ini  “  (d2  +  2dxr\  cos  d>'  +  z2  +  4/2)1/2. 


When  inserting  these  relations  into  Eq.  (11),  on  finds  that  the  current  J„  is  a  function  of  z,  and  <£' 
which  will  be  used  as  the  integration  variables.  The  vector  p  has  the  following  form: 

p  -  |(d,  +  r,  cos  4>')x  +  rx  sin  4>y  -  [ 4 f1  -  (d\  -  2dxrx  cos  <f>'  +  z,2)]z|/  |p|  (14) 


R  —  sin  8  cos  d>x  +  sin  0  sin  <I>^  +  cos  8r.  (15 

Therefore, 

p  ■  R  -  {(d(  +  rx  cos  <6')  sin  8  cos  <I>  +  rx  sin  d>'  sin  8  sin 

-  ~  [4/2  -  (d?  -  2d|Z,  cos$'  +  zjOcosSj/lpl.  (16 

Insert  this  relation  into  Eq.  (1)  and  perform  the  integration.  The  electric  field  has  the  following  form: 

E(S,  f  £7,0-,.  <f,')~expl-jk(p;-p  ■  R))r,dr,d<y.  (17 

n  Pi 


The  0  and  4>  components  are 


E%  -  cos  0(cos  <*>£■,  +  sin  4>Ey)  -  sin  <J»£, 
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and 

E#  -  -  sin  &EX  +  cos  (18b) 

where  Ex,  Ey ,  and  Ez  are  the  x,  y,  and  z  components  of  £( 8,  <J>)  in  Eq.  (17). 

COMPUTER  PROGRAM 

The  preceding  formulation  is  coded  into  a  FORTRAN  computer  program  which  has  been  tried  on 
the  Texas  Instruments  ASC  computer.  Appendix  B  lists  that  program.  The  FORTRAN  program  can 
compute  the  ratiation  pattern  of  a  reflector  which  is  only  a  portion  of  paraboloid.  This  portion  is  a  cir¬ 
cle,  has  a  radius  r{  and  center  at  d \  as  shown  in  Fig.  1 .  The  integration  is  then  computed  by  integrating 
the  <f>  angle  from  0  to  2 it  and  the  radius  rt  from  0  to  a  desired  value.  The  first  data  card  has  a  IS  fixed 
point  format  which  specifies  the  following  parameters: 

•  IEND  is  the  number  of  points  to  be  integrated  in  the  radius  r,  direction 

•  JEND  is  the  number  of  points  to  be  integrated  in  the  <£'  direction  (from  0  to  2 w  radians). 

•  NEL  is  the  number  of  feed  array  elements. 

•  NPT  is  the  number  of  feed  point  to  be  computed  at  a  given  <t>  angle. 

•  NPOL  specifies  the  feed  element  polarization,  if  NPOL  =  1,  the  excitation  field  is  polarized 
at  ^direction.  If  NPOL  =  0,  the  polarization  is  at  x-direction. 

The  second  data  card  has  an  8F10.6  floating  point  format. 

•  FO  is  the  focal  length. 

•  RSO  is  the  radius  of  the  complete  paraboloid. 

•  D1  is  the  center  of  the  circle  of  the  portion  of  the  paraboloid  to  be  used. 

•  A1  is  the  radius  of  the  circle  of  the  portion  of  the  paraboloid  to  be  used. 

•  FIC  is  the  <f>  angle  in  degrees. 

•  ARAG  is  the  range  of  ©  angel  to  be  plotted  in  degrees. 

All  of  the  preceding  data,  except  angles,  are  in  terms  of  wavelength. 

The  next  few  data  cards  specify  the  location  of  the  array  elements  in  x,  y,  and  z  directions.  These 
cards  are  all  in  F10.6  form;  each  contains  8  data.  Therefore,  the  number  of  cards  required  depends  on 
the  number  of  array  elements.  For  example,  if  the  array  has  24  elements,  three  cards  are  required  to 
specify  the  locations  of  each  coordinate.  Nine  cards  are  then  needed.  The  sequence  is  that  first  read  in 
the  x-component  locations,  then  the  y,  finally  the  z-locations.  The  output  of  this  program  is  a  plotted 
antenna  pattern  for  E%  and  £$  components  at  a  cut  of  <t>  angle  specified  by  FIC  and  plotted  a  range 
from  0  degrees  to  ARAG  degrees.  However,  if  ARAG  is  negative,  then  it  will  plot  from  -ARAG  to 
+  ARAG. 

A  function  subroutine  PATF  is  included  in  this  program.  This  function  specifies  the  array  ele¬ 
ment  pattern  in  its  own  coordinates  sin  0\  and  cos  d>'.  At  the  presented  time,  the  subroutine  function 
is  set  at  an  isotropic  element  pattern.  The  user  may  change  it  to  fit  his  requirement. 
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EXAMPLES 


Figures  2a  and  2b  show  the  radiation  pattern  of  a  paraboloid  dish  with  a  13.88-wavelength  radius 
and  a  25. 9- wavelength  focal  length.  A  single  source  located  at  the  focal  point  is  used  in  both  plots. 
Figure  2a  is  plotted  by  use  of  the  method  developed  by  Galindo-Israel  [21.  This  method  uses  a  Jaco¬ 
bian  series  to  expand  the  integral  of  Eq.  (1)  into  a  series  of  functions  in  terms  of  and  0.  This 
approach  should  be  efficient  in  computation  time.  However,  when  the  series  is  expanded  for  better 
convergence,  the  transformation  must  be  performed  at  the  vicinity  of  the  center  of  the  main  beam. 
Since  the  main  beam  position  is  a  function  of  the  sources’  location,  this  program  is  not  necessarily 
efficient  when  multiple  sources  are  used.  Figure  2b  is  plotted  by  use  of  the  method  described  in  the 
Mathematical  Formulation  section  of  this  report.  The  patterns  on  both  figures  are  exactly  the  same. 
The  central  processing  time  for  Fig.  2a  is  41.6  s  and  that  for  Fig.  2b  is  52.46  s.  Galindo-Israel’s  method 
is  about  30%  more  efficient. 


Fig.  2a  —  Radiation  pattern  of  a  paraboloid  reflector  computed 
by  the  series  expansion  method 


Fig.  2b  —  Radiation  pattern  of  a  paraboloid  reflector  computed 
by  the  exact  solution 
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Next,  we  tried  an  offset  case.  The  results  are  shown  in  Fig.  3.  The  required  central  process  times 
for  our  program  and  the  Galindo- Israel  program  are  about  the  same  as  that  in  Figs.  2a  and  2b.  Two 
cases  of  a  two-source  and  a  five-source  are  tried.  Their  patterns  are  shown  in  Fig.  4.  The  central  pro¬ 
cess  time  is  54.73  s  for  the  two-source  case  and  57.68  s  for  the  five-source  case.  There  is  no  way  to 
compare  with  Galindo-lsrael's  method  because  that  program  is  not  designed  for  multiple-feed  sources. 

From  these  examples,  the  increase  of  computer  time  for  multiple  sources  is  reasonable.  This  pro¬ 
gram  probably  is  efficient  in  this  sense. 


Fig.  3  —  Radiation  pattern  of  a  paraboloid 
reflector-off-set  feed 


Fig.  4a  —  Radiation  pattern  of  a  paraboloid 
reflector— two  multipie-source  off-set  feed 


CONCLUSIONS 

A  computer  program  which  computes  the  radiation  pattern  of  a  paraboloid  reflector  with  multiple 
feed  sources  is  presented.  From  some  computer  exampes,  the  program  is  efficient  for  a  large  number 
of  feed  sources. 
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Appendix  A 

DERIVATION  OF  EQUATION  (11) 


Pn  =  P'nzX  +  p'„yy  +  fal, 

+  unyy  +  umz. 
n  =•  nx  x  +  nyy  +  n2z, 

where 

P«  +  P'ny  +  P«  -  1. 

+  »ny  +  ui  -  1, 

nx  +  ny  +  n}  —  1 , 

and  veciors  fa ,  and  n  are  defined  in  Eq.  (5),  Eq.  (9),  and  Eq.  (10)  respectively. 

3,  “  Pn  X  «n  *  pi,. 

a  *  |  r  a  i  a  \ 

“  -  PnOn  '  «„)• 


Let 


Pn  '  ^  unxPnx  PiyiPiy  "1"  unxPnz  )  “  tt» 

«n  "  («w  -  <*P«)*  +  (Uny  -  «pi,y)y  +  (un  ~  «P*)^. 

e„  I  —  Vl  —  a2, 


*n 

7. 


r>*P«*e„ 

h  x  fax  uj-y/ 1  -  a2 

Pn  / «  «  \ 

-  (n  ■  «„)  - 


Vl  -  e2 


>/ 1  -  a  ‘ 


in  ■  fa) 


vl  — a2 


U-ww(n,p^  +  nzp«)  +  +  n„p^nz]x 


^ ^nxPny Pny^xP  nx  ^zPny)  PnzPny^z^y 


^nxPnz^x  PnyPny  Pz^xPnx  tyPny ) H 
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Appendix  B 

COMPUTER  PROGRAM  LIST 


CbN 


0  001 
0002 
0003 
000* 
0005 
0006 
0007 


0000 

0009 

0010 


oon 

0012 


0013 
001* 
0015 
0016 
0017 
0018 
0019 
OOcO 
0021 
0022 
0023 
002* 
002  5 
0026 
002  7 
00c  8 
002  9 
OOjO 
0031 
0032 
0033 
003* 


jOJRCc  LlbllNG  ASC  FAST  FORTRAN  COMPILER  KElEAS fc  FI=X05>»T 


STATcMcNT  CP  OPTIONS  =  CS.O.E.K.M.V)  OAT £  -  06/c*/82  CSN 


PROGRAM  REFPAT  0132 

;0MM0N/0/£c0-( b, 100)  0133 

COMPLEX  PF1,PP,PSUM,PC  013* 

DIMENSION  PFK3.500),PSUMC  3 . 50 0 >» P*U  3 * 5  00  >  . PP C  3)  »  J SPC  500  )  0135 

DIMENSION  ERAS<500),  aXC500),YU500>  0136 

3 1  ME  MS  ION  bI7AU200),COTAC(200),SIFI(200),COFI(200)  0137 

DIMENSION  tPLENC100),FJX(100),FJU100),Fj2C100)  0136 

C  THIS  PROGRAM  COMPUTES  THE  RADIATION  PATTERN  OF  OFF-SET  REFLECTOR  0139 

.  WITH  tXACT  SOLUTION.  RE : Cc C TOk  CAN  BE  PORTION  OF  TNt  REFLECTOR  01*0 

_  OF  A  PRABOLOID  AND  lT  CAN  HAVE  MORE  THEN  ONE  RADIATION  SOURCE  01*1 

CALL  RASTOP  01*2 

READ  IOO.IlND.JENO.NcLi  NPT.NPOL  01*3 

IOC  F0RMATC16I5)  01** 

C  IE NO,  NUMBER  OF  POINTS  INTEGRATED  IN  R  DIRECTION  01*5 

C  JENO,  NUMBER  OF  POINTS  INTEGRATED  IN  FI  DIRECTION  01*6 

.  IENO  AND  JENO  HUbT  8c  ODD  01*7 

MPT,  NJMBER  of  FIELO  POINTS  IN  THETA  DIRECTION  01*8 

C  NEL.  NUMBER  OF  ARRAY  ELEMENTS  01*9 

C  NP0L=1,  T  -POLAR 1 2 AT  ION , X=  0  >  X-P0LARI2AT ION  0150 

READ  101, F0.RS0.D1.A1.FIC. ARAB  0151 

101  F0RMATC8F10.6)  0152 

c  FO.  REFLECTOR  FOCAL  IENGTHCINMAVELENGTH)  0153 

C  RSO,  REFLECTOR  RADIUS  IN  NAVELEN&TH  015* 

C  31,  CENTER  OF  THc  CIRCLE  OF  PORTION  OF  THE  RcE LECTOR  UbED  0155 

C  Al,  RADIJS  Or  T He  JSZO  REFLECTOR  CIRCLE  0156 

C  FIC,  FIELO  PONT  AT  FIC  PLANE  015T 

C  ARA&  ANGLE  RANGE  °15b 

W  IF  ARAG.LT.O.  plot  -ARAG  TO  ARAG  0159 

PRINT  102,IEND,JENO,NPT,NEL,NPOL  0160 


102  FORM AT (10X, 1615) 

PRINT  103,F0,RS0t01,Al,FIC,  ARAG 

103  -ORNAT(10X,9-10.*) 

PI=3. 1*15926536 
ATR*PI/180. 

PI2*PI*2. 

r IC=I IC*RTA 
C OF IC-LObCFIC) 

SIFIC»SIN<FIC) 

XSL*8. 

»SL  =  5. 

IF(ARAG.GT.O.)GO  TO  j 
ARAGc«-ARAG*l. 

T  A- ARAGPATR 
AST-ARAG 

asf=-aras 

GO  TO  * 

3  ARA&2*  ARAG 
TA  =  0. 

AS  T  =  0. 

ASP* AR AG 
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AT 

SOJRCE  LIST  140  A SC  EAST  FORT 

CSN 

STATEMENT  CP  OPTIONS  = 

0035 

4 

TACINC=ARAG2/<NPr-13 

0036 

rAINC=rACI4C*ArR 

0047 

TAC=0» 

0038 

XSCC-XSC/ARA&2 

0039 

30  1  M~1,N>T 

0040 

PSUM(1»N>=CNPLX(0.«0.3 

004  1 

PSUMt2,N)=CNPLXC0.*0.) 

0042 

PSUMC3,N)=CNPLXtO. »0.) 

004  3 

XXCN3=7AC*XSCL 

0044 

COTACCN3-COSCTA3 

0045 

SITACCN)*SIN(TA3 

0046 

T  A  =  T  At T  AINC 

004  7 

1 

7 AC *7  AC* TAG INC 

0046 

TA-O. 

004  9 

TAlNC*PI//JtN0 

0050 

DO  2  J=1.JEN0 

0051 

COFKJ)  =  COSCTAJ 

0052 

SIFICJ>=$INCTA> 

005  3 

JN00=N0DCJ»2) 

0054 

JSPCJ)=JNOD42*(1-JNOO)*4 

0055 

I?<J.E4.1.3k. J.t3.J£ND)JSP(J)=l 

0056 

2 

l A=T  A*T  AINC 

005  7 

F02=F0*F0 

0058 

SlNC-Al/ IENO 

0059 

R1 =S INC 

C 

ENTER  OFF-SET  ARRAT  ELEMENT  LOCATIONS 

0060 

READ  101*CEL0CC1,I).1=1,NEL) 

0061 

READ  101 t (EL0C(2« I3>I-1>NEL) 

0062 

READ  101,(cL0CC3,n.I  =  l,NtO 

0063 

PRINT  103,C<£L0CCA,I).A*1.3),I=1.N£L) 

0064 

00  10  l=l,I£ND 

0065 

ROSO=Rl**2*Ol4»2 

0066 

R01*Rl*0l 

006  7 

30  11  N  =  1»  NPT 

0068 

PFICl«N>*CHPLXCO..O.) 

0069 

PFI(2,N>*CMPLXC0. .0.) 

0070 

11 

PFI<3,N)=CMPLXC0.,0.3 

0071 

30  20  J=t,JEND 

0072 

R*SQRTCRDSa*E.4R01*C3FlCJ)> 

0073 

RtOFI=C01*Rl*C0FKJ)) 

0074 

RSIFI*R1*SIFICJ) 

007  5 

R2=R4R 

0076 

F0R2*E0-R2/t4.»F0> 

0077 

12 

UN»S0RTCR2»4.*FO23 

0078 

PP(l)*CMPLXtO.,0.) 

0079 

PPC2J*CMPLXt0.,0.) 

0080 

PP<33«CHPLXL0.*0.3 

0081 

00  22  L* 1 >NEL 

0082 

PX*RCOFl-ELOCtliC> 

008  3 

PT=RSIFl-EL0tC2.L> 

0064 

Pt=-F0R2-EL0CC3,L) 

OOo  5 

»X2=PX**2 

0086 

PT2=PT**2 

0087 

P22*P2*»e 

0086 

EPLEN2sPX2*PT2*P22 

00B9 

tPU.'t<c>*SJRKt.PUN2J 

009  0 

IFCNP0L.CT.03C0  TO  21 

0091 

UM*S4RT<PT2*P22> 

009  2 

=  JN0R*UN4UM*tPLtNCO 

0093 

FJXCL)*<RSIFI4PT-2.*F0*P2>/FJM0R 

0094 

FJT(L)*~RC0F14PT /FJNOR 

0095 

FJ2CL)*-RC0FI*P2/FJN0R 

0096 

CO  TO  23 

0097 

21 

JM»SQRT<PX2*P22> 

0098 

FJN0R*UN*UM4EPLENCl) 

RELEASE  FTFXOS 
DAT  t  =  06/24/62 
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0099 
0100 
0101 
0102 
010  3 
0104 
0105 
0106 
0107 

0108 
0109 
0110 
0111 
0112 
0113 
0114 
0115 
0116 
0117 
0116 
0119 
012  0 
Old 
0122 
Old 
0124 

0126 

Old 

0127 

0128 

0129 

0130 

0131 


SOURCE  LISTING  AS-  EAST  FORTRAN  CONP  UEK 

iT  AT  ENEN  T  CP  OPTIONS  =  (6,U,E, K,N,Y3 

FJXtLl=-RSIFI*PX/FJNOR 
:JTCc)  =  CnC0FI*PX-2.»-0»P7)/p  JNOK 
F  JUL3  =  -RS1F1*PZ/FJN0R 
23  AUG=cPLtN<LJ*Plt 

PC*CNP  LX (LOSC  AUG 3  »SIN( AUG 3 ) 

PP(13=PP(13*- JX(L3*Pt 
PP(2)=PPC2)*F J»CL)*P- 
PPC33=PPC33*FJZCL3*PC 
22  CONTINUE 

-ORPUTc  »HASt  OUt  TO  DIFFERENT  FIElD  POINT 
RC  C  =  RCOF r»COF IC*RSIFI«SIFIC 
00  30  N* l , NP  T 

PHASE =P0R2  *L OT  AC <  N)-RCC*S I T ACCN3 
AUG-*HASt*F 12 

PC=CNPLX(C0SCAUG3. S INC  AUG)) 

PFICI»NJ=PF1C1»NJ*PPC13*JSPCJJ*PC 
PF IC2*N3=PFIC2,NJ*PPC2 )*JSPCJj*PC 
PF  K3,n)=P= IC3,N)*PP<3)*JS?(J)»PL 
30  CONTINUE 
2  U  CONTINUE 

IFCI.EO.  IEND)RSINP  =  R 
KSINP=jSPCjJ»R 
JO  16  L= li NP  r 

PSUNC1,L1*PSUNC1.L)*PFK1.LJ»RS1MP 
PSUNC2,  LJ*PSUN<2,LJ*PFU2,LJ*RSINP 
15  PSUNC3,L3=PSUH{3,l3*PF  It 3,l)*RsIRP 
10  R 1 =R l ♦ 3l NC 
C  FIND  PATTERN  FUNCTION 

PNOR*0. 

30  40  1*1,  NPI 

PC*  COTACCI)*CPSJRC1,I3*COFIC«PSUN<2, I )*S IF  1C >-PS JRT 3, 1)*SITAC(I3 
P ATC 1 , 1 3  =PC*C0NJGCPC3 
IF<PA7U,IJ.CT.PNOAJPNOA  =  PAT(1,JJ 
>C*-»SDNU,IJ*5I:lC*PSUH(2,Il«C0F  IC 
PAT  <  l , I ) = PC*L0NJG( PC ) 


RELEASt  F  T“X  05 
DATE  =  06/24/62 


12 


NRL  REPORT  8740 


»AT 

CSN 


0132 
0133 
0134 
0135 
0136 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
015  7 
0158 
0159 
0160 


SOURCE  LISTING 
STATENENT 


as;  fast  fortran  COAPILER  RELEASE  FTFX05 

CP  OPTIONS  =  (8  DATE  =  06/24/82 


IF  CP AT{2,I).GT.PNOk)PNOR=P ATC2.I) 

PRINT  105, CPSURU, I), 11  =  1,3) 

105  FORMATC10X,6E12.4) 

40  CONTINUE 

PRINT  104  » (P  AT  (1,I),I=1,NPT) 

PRINT  104, CPAT  C2, I),I=1,NPT) 

104  F0R«ATC//,E10X,8£12.4>) 

CALL  PLOTSCERAS.500,2.) 

CALL  ORIGIN<4.,0.) 

CALL  SQJARECO.,XSL,0.,TSL) 

CALL  NXAXIS(0.,0.,AST  ,  5., ASF  ,XSL. -.1 , . 1 , -1 ) 
CALL  NTAXISCO.,0.,-60. , 10. , 0. • VSL, 1 , . 1 , -1 > 

CALL  CENTER(XSL/2.,-.8,.3, 14NANGLE (BE GR EE S > . 0.,14) 
CALL  C£NT£R(-.5,TSL/2.,.3,15HA4D1ATI9N(DB>, 90. *13) 
00  50  Xs i, 2 
00  60  1  =  1  *  NP  I 

IF(PATCK,I).LE.0.)iO  TO  61 
OB=10.*ALOG10(PAT(K,I3/PNOR) 

TT CI )«C1. *08/60. >*TSL 
GO  TO  62 

6i  mn=-60. 

64  iF(mi).GT.rsLjrru3=rsL 
60  IF(YT(I).LT.0.m<I)  =  0. 

JJ  =  0 

IF(K.GT.1)JJ=10 

CALL  LINE(XX,YY,NPT,1,JJ,4) 

50  CONTINUE 

CALL  ENDPLT 
ENO 
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